Flat spin wave dispersion in a triangular antiferromagnet 
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The excitation spectrum of a S = 1/2 2D triangular quantum antiferromagnet is studied using 
1/S expansion. Due to the non-coUinearity of the classical ground state significant and non-trivial 
corrections to the spin wave spectrum appear already in the first order in 1/S in contrast to the 
square lattice antiferromagnet. The resulting magnon dispersion is almost fiat in a substantial 
portion of the Brillouin zone. Our results are in quantitative agreement with recent series expansion 
studies by Zheng, F jgirestad, Singh, McKenzie, and Coldea [Phys. Rev. Lett. 96, 057201 (2006) and 
|cond-mat/0608008|. 



Triangular antiferromagnets occupy a special niche in 
the studies of quantum magnetism. Ising antiferromag- 
net on triangular lattice has a finite zero-temperature 
entropy, which reflects an extensive degeneracy of the 
ground state manifold Classical Heisenberg model 
on a triangular lattice represents the textbook example 
of the full SU{2) symmetry breaking and noncoUinear 
spiral spin ordering in the ground state - the 120° spin 
structure. For a quantum S — 1/2 antiferromagnet on 
a triangular lattice, Anderson proposed back in 1973 the 
disordered resonating valence bond (RVB) ground state 
01 ■ This suggestion stimulated extensive research for 
over 25 years. An RVB ground state on a triangular 
lattice has been found recently [Sj], albeit for a quantum 
dimer model. It was also established, by large-N and 
gauge theory approaches, that a disordered ground state 
of a triangular antiferromagnet must possess unconfined 
massive spinon excitations 

On the experimental side, several novel materials with 
triangular structure have attracted substantial interest 
over the last few years. NaxCo02, in which Co atoms 
form a layered hexagonal structure [sl, was suggested to 
be in close proximity to a spin liquid [6| . Another poten- 
tial candidate for a spin liquid is k— {ET)2Cu2{CN)z at 
small pressures Ml . Finally, there is an intensive theoret- 
ical debate [Tfl [HI, on the structure of the ground 
state of a spatially anisotropic triangular S" = 1/2 anti- 
ferromagnet CS2CUCI4 [1]. 

The ideas about the disordered ground state of un- 
confined spinous, however, could not be immediately ap- 
plied to the most studied Heisenberg model of quantum 
S = 1/2 spins on a triangular lattice, as both perturba- 
tive 1/S and numerical calculations point that the clas- 
sical, 120° spin structure survives quantum fluctuations. 
Quantum fluctuations do reduce the average value of a 
sublattice magnetization to 50% of its classical value [l^ ■ 
This renormalization is generally comparable to that for 
a, S = 1/2 antiferromagnet on the square lattice fl^. For 
the latter, calculations to order l/S*^ for the spectrum do 
show that the overall scale of the spin-wave dispersion is 
renormalized by quantum fluctuations, but the disper- 



sion retains almost the same functional form as in the 
quasiclassical limit, and obviously is better described by 
magnons rather than by deconfined spinous. 



For a triangular antiferromagnet, Chubukov et al 15 1 
computed 1/S corrections to the two spin- wave velocities, 
and found that these corrections are quite small, even for 
S = 1/2. No 1/S' calculations of the full spin- wave dis- 
persion have been reported in the literature but, based 
on the results for the velocities, it was widely believed 
that the functional form of the dispersion for S = 1/2 
antiferromagnet should also be close to that in the qua- 
siclassical, large S hmit, i.e., that the spin- wave descrip- 
tion can be extended to S = 1/2. Recent series expan- 
sion study |13l] . however, uncovered remarkable changes 
in the functional form of the dispersion in the isotropic 
quantum S* = 1/2 antiferromagnet on triangular lattice 
compared to S" = 00 limit. In particular, the dispersion 
for the S" = 1/2 case possesses local minima ("rotons") 
at the mid-points of faces of the hexagonal Brillouin zone 
(BZ). The classical dispersion does not have such local 
minima. 

The authors of Ref. [3] conjectured that the qualita- 
tive changes between the actual dispersion for S = 1/2 
and the classical dispersion may imply that at energies 
comparable to the exchange integral J, the system is 
better described in terms of pairs of deconfined spinous 
rather than magnons (the latter in such description are 
bound states of spinous). In other words, they argued 
that the spinon description, valid for the disordered state, 
may adequately describe high-energy excitations of the 
ordered state. 

In this communication, we propose another explana- 
tion for the series expansion results, alternative to the 
one proposed in Ref. [l3|- We argue that regular 1/S 
corrections, extended to S* = 1/2, strongly modify the 
form of the magnon dispersion in a triangular antifer- 
romagnet, and the renormalized dispersion has "roton" 
minima at the faces of the BZ, in agreement with [l^. 
But we also found a more drastic effect - the renormalized 
dispersion turns out to be almost flat in a wide range of 
momenta. The flat renormalized dispersion has not been 
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FIG. 1: (Color online) BZ of a triangular antiferromagnet. 
Thin lines indicate the directions along which the "cuts" 
of dispersion are taken. Coordinates of the points are: 
P(27r/3, 0), A(7r, 0), Q(47r/3, 0), B(7r, tt/^S), C(27r/3, 27r/V3), 
F(0,27r/\/3), and E(0,7r/y3). ImE^®^ ^ inside the shaded 
star-shaped area. 



the decoupling of the four-boson term. These corrections 
do not affect the functional form of the dispersion [l^. 
The corrections to the dispersion then appear at l/S"^ or- 
der, from Tig term in the magnon Hamiltonian and from 
the second-order perturbation in 7^4 , both are very small 
numerically. In triangular antiferromagnets, Tis term is 
present due to the non-colUnearity of the classical spin 
configuration (120° structure) and, taken at the second 
order, gives rise to non-trivial corrections to the disper- 
sion already at the order 1/S. The expressions for Tis and 
7^4 have been obtained in [l5| and we refer to that work 
for the derivation. Decoupling the four-magnon term in 
a standard way and adding the result to the quadratic 
form we obtain, to order l/S, 



reported in [13| where the numerical data have been pre- 
sented only along special high-symmetry directions in the 
BZ. The subsequent, more detailed series expansion stud- 
ies 1J| , carried out in parallel with our research, did find 



the regions of flat dispersion. The analytical and series 
expansion results are in good qualitative and quantita- 
tive agreement, which indicates that the flat regions are 
likely present in the actual dispersion, and that a first or- 
der in 1/5* provides fairly accurate description of a trian- 
gular antiferromagnet even at high energies, comparable 
to the exchange J. 

The point of departure for the 1/S calculation is the 
Heisenberg antiferromagnet on a triangular lattice 



Ti = J ^ ^ S^Sj 



(1) 



<ij> 



Using Holstein-Primakoff transformation to bosons, as- 
suming the 120° spin structure in the ground state, and 
diagonalizing the__ quadratic form in bosons one can re- 
write Eq. H]) as 



n = 3JS 



L k 



(2) 



where dots stand for higher order terms in 1/S, and 



Ek = \ Al - Bl = v/(l-i^fc)(l + 2i/fe). 



(3) 



where 



Ak = 



Vk = 



Vk „ 3 



1 / I , n ^a; \/3ky 

- COS Kx + 2 COS — cos — - — 



(4) 



The BZ is presented in Fig. [T] The magnon energy Ek 
vanishes at the center of the zone, k — 0, where Vk — 
and at the corners of the BZ (e.g. at points Q and C), 
where i^k — —1/2. 

In a square-lattice antiferromagnet, Tis term is ab- 
sent, and 1/5 corrections to the dispersion come from 



n = 3JS 



Y^Ekclck + ^H: 



L k 



(5) 



where Ek represents dis per sion renormalized by quartic 
terms (coming from Hi jlSf l 



El^El{l + c,/S) + il-iyk)c2/S, 



(6) 



where ci = 1 - /q - /i/4 + 5/2/4, C2 = 3(/2 - /i)/4 and 
In = jf J2k ^k/^k- The cubic part reads 



1 

+ 3 



f 3 



\'i2N 



1/2 



^(cl4cfe'i>i(l,2,fc)(5 
1.2 



1+2-/C 



444*2(1, 2, k)5i+2+k + h.c. 



(7) 



Here 1 = ki, 2 = k2, and summation is over the BZ. The 
vertices $1,2(1, 2, k) = l'i,2(l, 2, k)/ ^/E[E^E^, where 

$1,2(1,2,3) = Pr/i^) (/i^Vf ± /i'Vi^') + ^2/i^^ X 

(^(1)^(3) ± ^(l)£3)^) ^ -^£3) ^^(1)^^(2) _ jWjmy^ 



are expressed in terms of /j_ — ■\/ . ± B^^ and 

\f%ky \ 



Vh — — sm — cos cos ■ 

3 2 I 2 2 



(9) 



The H3 term gives rise to a ^-dependent magnon self- 
energy to order 1/5 



1*1(1, 2, A;)|^ 



16SN\ ^ El 

l+2=k ^ 



E 

l+2=-k 



E2 — Ek 
*2(1,2,A:)|2 



iO 



El + E2 + Ek 



iO 



(10) 



Note that at this order $1,2 and Ek are expressed via 
bare (5 — 00) quantities. Collecting the corrections from 
three-magnon and four-magnon processes and restoring 
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the pre-factor 3JS (see Eq. [5]), we obtain for the full 
renormalized dispersion 



^,.e„(fc) = 3JS 



El 



(S) 



(11) 



Expanding (fTTj) to first order in l/S* we finally obtain 



Ek{l 



2S' 



(1 - ;^fc)c2 
2SEk 



(12) 

It was shown in jl5| that the renormalized dispersion 
preserves zeros at fc = and at the corners of the BZ, 
i.e., locations of Goldstone modes are not affected by 
1/S corrections. Below we compute the full spin- wave 
dispersion (fT^ and apply the results to S = 1/2. 

We used the Mathematica© software to calculate the 
self-energy E^, pop . Two-dimensional momentum inte- 
grals were regularized by replacing energy denominators 
with E + iS, where 5 = 10^'', and taking the real part of 
the resulting expression. The imaginary part was calcu- 
lated using the "bell approximation" for the delta func- 
tion: S{x) = ■\/i/7re^*^ . We found that t = lO'^ gives 
stable and consistent results. 

Our results are shown as three-dimensional plots in 
Fig El For comparison we also plotted the classical dis- 
persion SJSEk — 1.5JEk- For clarity, we only plot the 
dispersion over a quarter of the BZ, and set it to zero 
outside this quarter. Observe that the renormalized dis- 
persion vanishes at fc = and at the corners of the BZ 
hexagon. Spin wave velocity at fc = decreases in com- 
parison to the classical result, while the one at the corners 
of the BZ increases (see Fig 4(a) below) T^. 



We clearly see from Fig E] that the actual dispersion is 
rather different from l.bEk- The key difference is that the 
renormalized dispersion Eicn{k) has a plateau at around 
0.8 J over a wide range of momenta. This is most clearly 
seen in Fig 2(c) and |2(d)] The dispersion also possesses 
a roton-like minimum near the faces of the BZ, as is best 
seen along the cut FC, Fig |3(a) 



In Fig 2(e) and 2(f) we show the imaginary part of 



Eicn{k) from three- magnon processes. In distinction to a 
square-lattice antiferromagnet, the imaginary part of the 
dispersion in our case appears already at the order 1/S. 
As follows from PH)) . it is present when one-particle ex- 
citation (magnon) with momentum k can decay into the 
two-particle continuum, i.e., when -E^ + E^_^ = E^ for 
some q £ BZ. This condition defines star-shaped region, 



shown in light gray shading in Fig[T] (see also Fig 2(f) I 



inside it ImE^, is nonzero. While ImS^"^ ^ (and, 
hence, ImEronik) ^ 0) in most of the BZ, we found that 
Im_Eron(fc) does not exceed 0.26J, and is much smaller 
than Rei?ren(fc)- Nonetheless, a finite imaginary part is 
important as it responsible for dampening excitations at 
wavevectors where variations of Rei?ron(fc) in Fig 2(c) 



and 2(d) are maximal. 

Observe that the maxima of Im£'ren(^) in Fig 2(e) 
occur exactly where the variation of Re-Eicn(fc) is the 
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FIG. 2: (Color online) Classical and renormalized spin wave 
dispersions. Left panels - three-dimensional plots, right pan- 
els - color density plots. The dispersions are shown in the 
OQCF quadrant of the BZ and are set to zero outside this 
region to ease the viewing. The classical dispersion, figures 
(a,b), is 1.5i5fe. The real part of renormalized dispersion 
Ercn{k) (Eq. WI\i is shown in figures (c) and (d). The flat 
dispersion is in the blue (dark) region in (d). The imaginary 
part of the renormalized dispersion is shown in figures (e, f). 
Observe that it is numerically small throughout the BZ, and 
is zero outside the "star" region, see also Fig[T] 



strongest, see also Figs [3][S] below. In other words, the 
sharp features in the 1/S dispersion are not artefacts of 
numerical calculations, but are real features of the excita- 
tion spectrum at this order. Observe also that Imi?i.en(^) 
vanishes at momenta where ReSron {k) stays almost con- 
stant. This implies that nearly immobile magnons have 
infinite lifetime, i.e., are true excited states of the system. 

To better display the new features and for comparison 
with series expansion studies, we show in Figs[5][5]the dis- 
persion E^on{k) along six different cuts through the BZ 
(black lines). The directions of particular cuts are shown 
in Fig[T] In each of the plots, we also show the classical 
dispersion (dashed red lines) . The flat regions are clearly 
seen in the cuts along FC, OF, PC, and EB directions, 
and the roton minimum is seen in the cut along FC di- 
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FIG. 3: (Color online) Classical and renormalized dispersions 
(ordinate) along (a) FC and (b) OF directions (abscissa). Red 
dashed line - classical dispersion. Black solid line - the real 
part of the renormalized dispersion Eren{k), Eq. (|12|l . Blue 
dotted line - the imaginary part of E^euik). The black dots 
are series expansion data from [ij. No fitting parameter is 
involved in the comparison with the data. 
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FIG. 4: (Color online) Same as in Fig. |3] but along (a) OQ 
and (b) PC directions. 



rection. Note that in accordance with Fig [1] the local 



minimum near point F (Fig 3(a) I corresponds to a truly 
stable excitation, Im£'ren(^) = there. The maximum of 
the renormalized dispersion is at 1.2J which is substan- 
tially smaller than 1.6J for the classical dispersion. 

We now compare in some detail our results to series 
expansion studies of Zheng et al [l^ 14 1. First, series 



expansion studies have found that the dispersion has a 
bandwidth of about J (see series data in Fig |4(b)[ ), which 
is substantially smaller than 1.6 J for a classical disper- 
sion. This is in agreement with our results (see Figs.[3][5]). 
Note that this difference is much larger than one might 
expect by comparing the spin-wave velocities which are 
renormalized only by about 10% [15|. Second, series ex- 
pansion studies found "roton" minima near points F and 
B in the BZ, Fig[T] Our cuts along FC and OF also show 
a minumum near point F (see Fig [3]). Third, recent series 
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FIG. 5: (Color online) Same as in Fig. [3] but along (a) AB 
and (b) EB directions. 



expansion studies fT3| found the regions of flat dispersion 
at 0.8 J along FC, PC, and EB directions. Our results for 
E^cn{k) also show the regions of flat dispersion at around 
0.8 — 0.9 J. The flat regions are clearly seen on the den- 
sity plot in Fig 2(d)[ as well as in Figsjg |4(b)] and m 
Fig |5(b)| (near points E and B). 

The near-flat regions of the magnon dispersion have a 
large density of states and can be probed by Raman scat- 
tering. The two-magnon Raman intensity in a square- 
lattice antiferromagnet is peaked slightly below twice the 
frequency at which the magnon density of states diverges. 
In our case, the density of states is peaked at around 0.8 J, 
and we expect that the two-magnon Raman intensity will 
be peaked somewhat below 1.6J. 

To conclude, in this paper we used 1/5" expansion, ex- 
tended it to S — 1/2, and obtained the renormalized 
magnon dispersion for a Heisenberg antiferromagnet on 
triangular lattice. We found that the renormalized dis- 
persion is qualitatively different from the classical one - 
it is almost flat over a wide range of momenta in the BZ, 
and has roton-like minuma near the mid-points of faces 
of the BZ. These results are in full agreement with recent 
series expansion studies. 

We are thankful to R.R.P. Singh for valuable discus- 
sion and to R.R.P. Singh and W. Zheng for sending us 
series expansion results prior to publication. OAS ac- 
knowledges the donors of the American Chemical Soci- 
ety Petroleum Research Fund for support (PRF43219- 
ACIO). Ave and AGA are supported by NSF via grants 
No. DMR-0240238 and DMR-0348358, respectively. The 
authors thank the Theory Institute at Brookhaven Na- 
tional Laboratory where this work has initiated. 

JVote added: after submission of our manuscript we 
learned of a closely related study of magnon decays in 
triangular antiferromagnets (l7| . 
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